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Dynamical characteristics of Rydberg electrons released by a weak electric field 


Elias Diesen, 1 Ulf Saalmann, 1 Martin Richter, 2 Maksim Kunitski, 2 Reinhard Dorner, 2 and Jan M. Rost 1,3 

1 Max-Planck-Institut fur Physik komplexer Systeme, Nothnitzer Strafle 38, 01187 Dresden, Germany 
2 Institut fur Kernphysik, Goethe-Universitat, Max-von-Laue-Strafle 1, 60438 Frankfurt am Main, Germany 
3 PULSE Institute, Stanford University and SLAC National Accelerator Laboratory, 

2575 Sand Hill Road, Menlo Park, California 94025, USA 

The dynamics of ultra-slow electrons in the combined potential of an ionic core and a static electric 
field is discussed. With state-of-the-art detection it is possible to create such electrons through strong 
intense-field photo-absorption and to detect them via high-resolution time-of-flight spectroscopy 
despite their very low kinetic energy. The characteristic feature of their momentum spectrum, 
which emerges at the same position for different laser orientations, is derived and could be revealed 
experimentally with an energy resolution of the order of 1 meV. 
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Extracting excited electrons from an atom or molecule, 
even a plasma with a constant electric field is an 
established technique which can reveal the minimal 
(Rydberg)-excitation by tracking the electron yield as a 
function of applied field strength JT] . In the time domain, 
very small extraction fields of about F = 1... 10 V/cm 
(with 5.142 V/cm = 10 -9 atomic units) and high, pulsed 
Rydberg excitation lead to intricate electron dynamics 
despite the fact that a hydrogenic problem in an electric 
field is separable (e.g., in semi-parabolic coordinates). 
Clearly, electrons which are capable to escape under such 
conditions must be highly excited and this is achieved in 
the experiment by a preceding excitation with a short 
intense laser pulse [2]. Interestingly, the details of this 
excitation are irrelevant for the dynamics described here. 

The relevant questions of this problem are: Can one 
describe the dynamics classically or semi-classically? We 
will demonstrate, in comparison to experiment, that a 
classical description is very accurate and we will explain 
why this is the case. A second question regards the mo¬ 
mentum spectrum of these electrons: Is it structured or 
smooth? We will demonstrate that it is characterized 
by a main peak, offset from zero by a well defined mo¬ 
mentum, in very good agreement with experiment. The 
dependence of the peak position on the field strength F 
is deduced analytically and compared to experimental 
results without any free parameters. 

The kind of Stark dynamics discussed here with a fo¬ 
cus on the differential momentum distribution of ionized 
electrons has not been investigated before, since ioniza¬ 
tion of Rydberg states in static or pulsed weak electric 
fields has mainly served the purpose to extract details 
of specific quantum states ranging up to principal quan¬ 
tum number n ~ 100 by either tunneling or over-barrier 
ionization (ZEKE spectroscopy) [l][3j[4]. Individual clas¬ 
sical trajectories and their contributions to the electron 
yield, on the other hand, have been investigated both in 
the context of astrophysics, where the same Stark Hamil¬ 
tonian arises from a combined gravitational and constant 
driving field [5], and in an atomic setting. For the latter, 
the inclusion of phases to account for interference effects 
of different pathways to a detector m even lead beyond 


a classical treatment. 

The results presented here are also relevant for all 
experiments using electric-field extraction techniques, 
in particular the COLTRIMS [9] technique which is 
nowadays used routinely in many reaction microscopes 
(REMI) world wide. At very low energies, the Stark dy¬ 
namics described here will modify the results expected in 
such setups. Moreover, it adds another interesting fea¬ 
ture to the low-energy kinetic-energy spectrum of elec¬ 
trons ionized by an intense laser pulse, a process which 
has received a lot of attention [MIS] since the discovery 
of the so called low-energy structure mm- 

To understand the phenomena under discussion, it is 
sufficient to consider a hydrogenic problem whose dynam¬ 
ics is governed by the Hamiltonian (atomic units are used 
unless stated otherwise and bold letters describe vectors) 

Pi = H + r • f (t) cos cut (1) 

with the Stark Hamiltonian 

H = —-— — Fz, r = \/x 2 +y 2 +z 2 (2) 

2 r 

where the ionic charge Q is taken for the rest of the pa¬ 
per to be Q = 1 (but could assume also an effective non¬ 
integer value) and a constant electric field of strength F 
points in the negative z-direction. The second term in 0 
describes the interaction with the laser pulse of (strong 
but otherwise arbitrary) envelope f(t), linearly polarized 
along or perpendicular to the direction of the extraction 
field along the z-direction in H. It is this combination 
of fields which leads to the “zero-energy” structure ob¬ 
served in the strong field photo ionization spectrum with 
REMIs [IS]. However, the role of the strong laser pulse 
is merely to populate the Rydberg states which serve as 
an initial condition for the unusual Stark dynamics we 
are going to discuss. Any other means to create the Ry¬ 
dberg population would lead to the same result if the 
Rydberg atoms are placed in a tiny electric field of the 
order of 10 V/cm: The photo-electron spectrum exhibits 
a peak not at zero energy but slightly upwards at a field- 
corrected photo-electron momentum (defined in Eq. (|3| 
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FIG. 1: a) 2D momentum map of photo-electrons from Ar 
atoms ionized with a laser pulse f(t) = e z f(t/T ) cos(cjt) with 
T = 64 fs and uj = 0.3875 eV (A = 3.2 /im) and an extraction 
field F = 1.5V/cm along the z-axis. 

Lower panels: Spectra as obtained by integrating the 2D dis¬ 
tributions over b) p_L = 0 ... 0.35, c) p_L = 0.02 ... 0.06, and 
d) p^ = 0.0.. . 0.02, respectively. Calculations (red lines) are 
compared to experimental data (blue, [18]). Note that these 
spectra are symmetrized with respect to IIy. 

below) of n* = — 0.6F 1 / 4 , with the asterisk indicating 
the peak value. 

Note that we have assumed i z = 0 (equivalent to the 
magnetic quantum number m = 0) and thus ignore the 
dynamics in the cyclic azimuthal angle. This reduction is 
possible as the laser-induced dynamics occurs on a length 
scale (given by the quiver amplitude, typically a few nm), 
which is much smaller than the barrier distance (typically 
fi m) of the Stark Hamiltonian and confirmed by the 
excellent agreement with full 3D simulations. Since the 
momentum of the electron increases under H as p z oc 
or p z oc t, with distance z from the nucleus or time t, 
respectively, its value depends on where it is measured. 
To avoid such a dependence, the experimental time-of- 
flight t Q f is used to construct an effective momentum 

n = p z (t 0 f) — Ft of , (3) 

which eliminates the effect of the constant electric field 
asymptotically. In other words, n does not depend on 
the position of the detector as long as it is placed suffi¬ 
ciently far beyond the barrier. This is a convenient pro¬ 
cedure since experimentally the time of flight t Q f and the 
detector distance Zd are known. Consequently, one ob¬ 
tains n = Zd/t 0 f — Ft 0 f/2, which is equivalent to ©• In 
quantum mechanics, the momentum variable n can be 
elegantly obtained directly through the interaction rep¬ 
resentation with H t = UHU^ where U = exp(—lFzt). 
Similarly, in classical mechanics, a canonical transforma¬ 
tion with the generating function — (p z +t)z leads to 


new canonical variables p z = n and z = z. 

We note that in all experimental studies using 
COLTRIMS reaction microscopes n in Eq. is referred 
to simply as the momentum component p z of the elec¬ 
tron parallel to the field direction of the spectrometer [9!. 
The underlying tacit assumption is that the motion of the 
electron can be split into an interaction region, where the 
extraction field can be neglected, and an asymptotic re¬ 
gion, into which the electron is launched with an initial 
momentum p z . It is this separation which breaks down 
for the situation discussed here. 

In Fig.[l] we compare classical simulation according 
to Fi given in ([!]) with experimental results [18] ob¬ 
tained for argon atoms in a strong few-cycle pulse with 
A = 3200 nm. Those experiments were the first to re¬ 
port a “zero-peak”, i. e. an extremely narrow contribu¬ 
tion of near-zero momentum electrons m, subsequently 
confirmed and attributed to Rydberg electrons [T9] . We 
use standard tunnel-ionization probabilities mm and 
propagate electrons according to Newton equations of 
motion for about 10 7 trajectories with the Hamiltonian 
in the attractive Coulomb potential, the driving 
laser pulse f (t) cos cot and the extraction field F. Hence, 
the calculations comprise the formation of Rydberg elec¬ 
trons [2] as well the subsequent Stark dynamics. The 
experimental results [181 shown in Fig.[l]D-d are sym¬ 
metrized in n with respect to n = 0 and show a double 
peak very close to n = 0. As one can see, the likewise 
symmetrized classical spectra agree extremely well with 
the experiment. 

However, a careful inspection of the contour plot in 
Fig# reveals that the Rydberg peak is asymmetric. In 
fact the peak position n* has a negative value, pointing 
into the opposite direction of the extraction field F, i. e., 
parallel to the actual electric field. This is confirmed by 
our experiment shown in Fig. [2] where for a much shorter 
wavelength of the driving laser but with the same extrac¬ 
tion field a similar peak in the momentum distribution is 
observed, shifted from zero. Switching to a setup, where 
the laser polarization is along the y- axis, i. e., perpen¬ 
dicular to the extraction field along z, finally confirms 





FIG. 2: Measured longitudinal momentum distribution for N 2 
molecules ionized with an 780 nm pulse with duration of 45 fs. 
As in Fig.[l] the polarization of the laser pulse is parallel to 
the extraction field of F — 2.6V/cm. 














3 


that the peak and its negative offset in n is exclusively 
due to the extraction field F, see Fig-U Both the calcu¬ 
lation and the experiment reveal the peak close to zero 
and again with a negative shift in the static field direc¬ 
tion, here n_L- We also note the qualitative agreement 
between theory (for argon atoms) and experiment (for 
N 2 molecules) which demonstrates that the laser pulse 
only prepares the initial Rydberg distribution for the en¬ 
suing Stark-held dynamics. It takes place far away from 
the (ionic) core whose exact nature hardly matters. 

Backed by these results we are now in a position to con¬ 
centrate on the Hamiltonian of © classically. In a first 
step we exploit its scaling properties with scaled phase- 
space variables 

{x, y, z, r} = F 1/2 {x, y, z, r} (4a) 

{Px,Py,Pz} = F~ 1/4 { p x ,p y ,pz} (4b) 


and scaled time t = F 3 / 4 1, which eliminate the depen¬ 
dence on F in Q with the new Hamiltonian 



( 5 ) 


This has two advantages: Firstly, the (experimental) de¬ 
pendence of certain observables on the held strength F 
can be immediately predicted, e. g., p z = F 1 / 4 ^, which 
applies to all momenta including the effective n = F 1 / 4 ^ 
with 7 r=p z — t. Secondly, the numerical calculation is 
largely facilitated since it is reduced to the problem of 
a single effective held strength (F=l) with a reasonable 



FIG. 3: Momentum distribution of photoelectrons with 
|p x |<0.01, whereby the extraction held (2.6V/cm) is perpen¬ 
dicular to the polarisation of the laser pulse (780 nm), i. e. 
f_Le^, which is in contrast to Figs, [l] and [2] Upper panel: 
calculation for Ar. Lower panel: measurement for N 2 . Each 
panel shows the 2D distribution and the one integrated over 
the parallel momentum in the interval |p y | < 0.005. 
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FIG. 4: Classical momentum spectra generated from a uni¬ 
form population of energies E s= —2... 0 and an isotropic 
distribution of cos ft n it • a) Contour plot as a function of the 
scaled hnal momenta 7r and p Q . The dashed lines show the 
hnal momenta for 4 energies indicated in the graph and initial 
angles in the range ft n it = 0 ... arccos(F 2 /2—1). The arrows 
indicate the location of the zeros Pg(7r n )—0 for E — 0 accord¬ 
ing to 0- b) Spectrum as a function of 7r after integration 
over the radial momentum \p Q \ < 0.1. 


(effectively an atomic-scale) range for the phase-space 
variables. 

From the scaling we immediately deduce that the off¬ 
set of the peak position n* from zero should scale with 
F 1 / 4 . Indeed, the experiment confirms this dependence 
as shown in Fig. [5] While the shift of the peak position 
with increasing extraction voltage or extraction field, re¬ 
spectively, might have been expected, the very existence 
of a peak is not obvious on a first glance since one would 
expect that the (almost homogeneously) distributed Ryd¬ 
berg population is set free by the extraction field in a cer¬ 
tain interval AF ... 0, with the lower bound A E given by 
the barrier formed by the static field. This should lead 
to a finite value of the photo-electron spectrum at zero 
energy but not to a peak, even shifted from zero in the 
observable n*. 

To gain insight into the origin of the peak and its po¬ 
sition 7 r* we take a closer look at the classical dynamics 
of the Stark Hamiltonian ([5|. To this end we propagate 
trajectories starting at the origin z = g = 0, with the 
“directed” distance g = ±^/x 2 -\-y 2 from the z- axis. The 
sign is positive (negative) if the electron has passed the 
z-axis an even (odd) number of times [21]. The initial 
distribution is uniform in energy F and angle cos#i n i t 
in accordance with the uniform occupation of Rydberg 
states. It turns out (as most easily seen in semi-parabolic 
coordinates, discussed below and in the supplement [22]), 
that for E < 0 initial angles cos#i n i t < F 2 /2—1 lead to 
bound trajectories, despite the fact that they have an 
energy above the barrier Fb = —2. These are the well- 













4 


known bound states in the continuum pQ, which are a 
consequence of the separability of the hydrogen problem 
in an electric field [5, 23] and account for 1/3 of all tra¬ 
jectories in the range E = —2 ... 0. The other 2/3 of 
the trajectories will escape. The abundance of the re¬ 
sulting final momenta 7r and p Q for those free trajectories 
is shown as contour plot in Fig.[4}i. The final momen¬ 
tum p Q = Ty/pcc 2 + p y 2 perpendicular to the extraction 
field oscillates as a function of 7r, with decreasing ampli¬ 
tude and period for each fixed energy E (dashed lines). 
Amplitudes also decrease at a fixed momentum p Q for de¬ 
creasing energy E. This topology gives rise to a caustic 
of final momenta around {tt *,p*} = {—0.6,0}. 

The peculiar structure of the momentum map in 
Fig.|4}i becomes understandable when taking advantage 
of the separability of the Stark Hamiltonian in semi¬ 
parabolic coordinates [5l [23] {u, v} = {y/r+z, Ey/r—z}. 
Here the sign for v corresponds to the sign of g. In these 
coordinates H from © assumes the form [22 

p 2 /2 - Eu 2 - u 4 /2 = l + /3, (6a) 

p 2 / 2 — Ev 2 + v 4 /2 = 1-/3, (6b) 


where (3 is the separation constant due to the addi¬ 
tional dynamical symmetry. This separation constant 
is fixed by the initial angle (in real space) according to 
f3 = cos #i n it as shown in the supplement [22]. The sepa¬ 
ration © is connected with a new time r which is related 
to the real time t via d t = (u 2 +v 2 ) dr. 

While the separability of the dynamics in semi- 
parabolic coordinates is of no significant help to produce 
the momentum map in Fig.[4}i, it does allow one to un¬ 
derstand the origin of the oscillating final electron mo¬ 
mentum 7 r(p g ). First we note [22 that asymptotically 
for large z or u the cartesian momenta are mapped as 
{PziPq} {u,v}. As can be seen from (p|, the dy¬ 
namics is that of a bound (v) and inverted (u) quartic 
oscillator at an energy 1 =f jd. The partition of the ef¬ 
fective energy between the two modes is controlled by 
/3. For f3 = 1 or cos#i n i t = 0 all energy is in u leading 
to the quickest possible escape with the largest value i r 
for u oc while no energy is left for the bound oscil¬ 
lator which stays at the fixed point v = 0. Decreasing 
/3, keeping E fixed, reduces the effective energy in u and 
provides in turn some energy for oscillations in v. Once 
there is enough energy for the v oscillator to complete 
half an oscillation, while the electron escapes, it will end 
up again at i?(t—^Tqo) = 0. Hereby is the escape time 
in ^-direction, which is finite due to the quartic term in 
(6a). Thus we can define a series of 7r n , for which asymp¬ 
totically v = 0 and therefore p 6 = 0, with n = 0,1,... 
the number of half oscillations in v. That these values 7r n 
are finite despite the ever increasing momentum due to 
the acceleration from the electric field is the consequence 
of the definition of tt = p z — £, cf. ©, which compensates 
exactly for this acceleration. 

Analytical quadratures for the phase space variables 
of the quartic oscillators in form of combinations of el¬ 
liptic functions exist [5], but do not yield simple expres¬ 



F [V/cm] 

FIG. 5: Shifting of peak position n* with the extraction field 
strength F for photo-ionization of N 2 , cf. Fig.[2] Note that the 
pre-factor in II* = — 0.6F 1 / 4 has been obtained theoretically 
from Fig.[4]and is not a fitting factor. 


sions describing the prominent features of Fig. [4] Yet, 
for the highest energy E = 0 these integrals simplify and 
therefore it is possible to determine the zeros of p e ( tt) 
analytically, 


, _ 2(1-2n 2 ) aAF(3/4) 

k=° “ (1 + 4n 4 )!/ 4 T(l/4) : 


( 7 ) 


as derived in detail in the supplement [22] . The arrows 
in Fig. ©> indicate the positions of the 7r n according to 
0 in good agreement with the numerical results. 

While the 7r n do not give the exact clustering positions 
they follow the same qualitative evolution and provide 
analytical insight for the latter: The second clustering 
near 772 is much weaker because it takes place at smaller 
values of (3 as compared to 7Ti. For the zero crossings at 
E = 0 the corresponding values read /3 n = (1 — n 4 )/(l + 
n 4 ). Since for escape, i. e. over-barrier motion, E 2 /2 < 
1 + (3 the phase-space volume contributing to the caustics 
decreases with the order n and the corresponding peak 
in the spectrum gets weaker. 

The excellent agreement of our classical calculations 
with our experimental results raises the question why 
quantum effects do not play a role, in particular given the 
important role of the potential barrier for the present dy¬ 
namics. The reason is the weakness of the field which ren¬ 
ders any action like quantity very large as can be directly 
inferred from the scaling (4a), S = SF -1 / 4 . This ap¬ 


plies to the imaginary action of a tunnelling path, so that 
tunnelling and diffractive phenomena are comparatively 
weak. The same argument holds true also for the real 
actions appearing as phases. Since the action differences 
scale with F -1 / 4 as well and are large for small fields, in¬ 
terference between different trajectories is quenched leav¬ 
ing the classical limit as an excellent approximation. This 
is in contrast to the situation for strong fields where in¬ 
terference patterns are observed [24, 21], a t fields several 
orders of magnitude stronger than discussed here. 

In summary we have shown that the Stark dynamics of 
escaping Rydberg electrons leads to a peak in the spec- 
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trum of the effective momentum n or the time of flight 
as the experimental observable, changing its position 
with the field strength according to U*(F) = — 0 .6F 1 / 4 . 
This explains the occurrence of the “zero-energy struc¬ 
ture” |T8]. The universality of the phenomenon is un¬ 
derlined through the very good agreement with experi¬ 
ments where the initial Rydberg population was created 
by laser pulses of different wavelengths and for differ¬ 
ent targets (atoms and molecules). This peak is not only 


an unexpected and intricate consequence of the Stark dy¬ 
namics, it also has direct implications for the experiment: 
The peak can be used to gauge the momentum scale if 
the field is known or vice versa, to determine the field 
strength with known momentum scale. 
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